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Abstract 

In the subcritical interval of the Reynolds number 4320 < R < 
R c = 5772, the Navier-Stokes equations of the two-dimensional plane 
Poiseuille flow are approximated by a 22-dimensional Galerkin represen¬ 
tation formed from eigenfunctions of the Orr-Sommerfeld equation. The 
resulting dynamical system is brought into a generalized normal form 
which is characterized by a disposable parameter controlling the magni¬ 
tude of denominators of the normal form transformation. As rigorously 
proved, the generalized normal form decouples into a low-dimensional 
dominant and a slaved subsystem. From the dominant system the crit¬ 
ical amplitude is calculated as a function of the Reynolds number. As 
compared with the Landau method, which works down to R = 5300, 
the phase velocity of the critical mode agrees within 1%; the critical 
amplitude is reproduced similarly well except close to the critical point, 
where the maximal error is about 16%. We also examine boundary 
conditions which partly differ from the usual ones. 


1 Introduction 

In local bifurcation theory the center manifold [||, Landau [j2j. |3j, and normal 
form 0, i i methods have been invaluable tools. As a common feature, 
a dynamical system is approximated by series expansions around the fixed 
point in phase space, and in the case of the first two methods also around 
the critical point in parameter space. In many cases rigorous theorems are 
available to ensure the correct behaviour of the system in a sufficiently small 
neighbourhood of the critical point. 

In this article we present a generalized normal form scheme with the aim to 
explore the subcritical , nonlinear regime in a finite neighbourhood of the fixed 
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point. The method is characterized by a disposable parameter e which con¬ 
trols the magnitude of the denominators of the normal form transformation. 
The case e = 0 corresponds to the most usual normal form where resonance 
denominators are avoided. With e > 0 we additionally discard quasi-resonant 
denominators with absolute magnitude < e. While increasing e enlarges the 
definition domain of the transformation, the normal form gets more and more 
complicated by the occurrence of further nonlinear terms. Nevertheless, as 
rigorously proved under rather general conditions on e, the normal form de¬ 
couples into a low-dimensional dominant and a slaved |7j subsystem in the 
sense that stability of the dominant system induces stability of the remaining 
degrees of freedom. The theorem proved generalizes previous work . 

We apply the method to the well studied subcritical stability analysis of 
the plane Poiseuille flow, see || |TT]. [ll], [12], |13[ [L4|, [15], [l6| . The Navier-Stokes 


equations are approximated by a Galerkin representation in the basis of 22 
eigenfunctions of the linearized part, i.e. the Orr-Sommerfeld equation. The 
Galerkin system, with diagonal linear part, is then subject to a normal form 
transformation with the resonance condition generalized through the interval 
parameter e. The latter is chosen so as to obtain a three-dimensional dominant 
system for the normal form in the Reynolds number interval 4320 < R < R c . 
The dominant system exhibits a Hopf-bifurcation at a critical amplitude r c = 
r c (R). The normal form series expansion is carried out up to cubic monomials. 


As compared with previous investigations of the plane Poiseuille flow JTl 


II, [T(| we examine the effect of slightly different boundary conditions: for the 


mode with zero wave number q in downstream direction we adopt the same 
boundary conditions as for the other modes. Because the different modes are 
mixed by the nonlinear interaction, we think one should stay in the same 
function space for consistency. 

To establish the Orr-Sommerfeld eigenfunctions of the linear problem we 
adopt the usual scheme of Chebyshev polynomial expansions |17|. We use 50 
even and 50 odd polynomials. Eventually the critical curve r c (R) obtained 
from the normal form is transformed back into the original phase space and 
compared with the Landau method [|T2| . 

The definition domain of the method is principally determined by the con¬ 
dition that the Jacobian J of the normal form transformation differs from zero 
in an open neighbourhood of the fixed point, where J = 1. Because of our 
cut-off at cubic terms the corresponding J does not allow for a reliable ex¬ 
trapolation. In the given case we check the definition domain in a pragmatic 
way by integrating the Galerkin system directly for starting points around the 
critical curve found by the normal form method. 

We set a — a c — 1.02 always, which is the critical value of the basic wave 
number in the downstream direction at the critical Reynolds number R = R c . 

The article is organized as follows: In the following section the two-dimensional 
Navier-Stokes equations are expressed in terms of the scalar stream function 
T(:r, z). The numerical methods to obtain the eigenfunctions of the linearized 
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stationary part are briefly sketched and the new boundary conditions are 
stated. In section 3 the Galerkin representation is derived. In the next section 
our generalized normal form scheme is introduced and a decoupling theorem 
is formulated which is proved in Appendix B. The fifth section contains the 
application of the normal form scheme to the nonlinear stability analysis of the 
plane Poiseuillc flow. In particular, the detailed normal form with numerical 
coefficients is given at the lowest Reynolds number considered R = 4320. In 
section 6 the quantitative results are presented and compared with the results 
of the Landau method 


12fl . This is followed by a conclusion section. 


2 Basic scheme 


The two-dimensional velocity field (u, w ) of the incompressible plane Poiseuillc 
flow is most conveniently described in terms of the stream function T as 


u = 


<9T 

dz 


w = — 


<9T 

dx 


( 1 ) 


where u denotes the streamwise velocity component and w the component 
normal to the boundaries. The Navier-Stokes equation then reads |L8 ; 0 


d T <9T <9V 2 T 

—v 2 t h- 

dt dz dx 


3* dV 2 * = i v% 

dx dz R 


( 2 ) 


The stream function T is decomposed as T = T 6 + where the basic flow 
T 6 = z — (l/3)2; 3 and the disturbance field is written as a superposition of 
plane waves 

OO 

^ d (x,z,t) = E * g (z, t) exp(iqax) (3) 

q =—oo 

with basic wave number a. All magnitudes are dimensionless: the coordinates 
are measured in terms of the channel half-width h, (u, w ) by the maximal 
unperturbed velocity Uq. The Reynolds number R = U^h/v where v denotes 
the kinematic viscosity. With the ansatz (|3|), equation (pi) takes the following 
well-known form ||: 


d 


L^ q {z,t)-R-(D 2 


q 2 a 2 )ty q (z,t) = 


where 

with D 


OO 

-iRa J2 [(<? - q')y q -q'(z,t)(D 3 - q' 2 a 2 D)^ ql (z,t) 

q '=—oo 

- q'Dy q _ q ,(z,t)(D 2 - q ,2 a 2 )* q ,(z,t)\ , 

L = ( D 2 — q 2 a 2 ) 2 — iaqR [2 + U{z){D 2 — q 2 a 2 )} , 
d/dz and U(z) = 1 — z 2 . 


(4) 

(5) 
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The linearized, stationary part of (J|) with \P 9 (z, t) = <3>g(;r) exp(A g f) leads 
to the Orr-Sommerfeld equation 

L%{z) = \ q R(D 2 - q 2 a 2 )$ q (z) . (6) 

We adopt the boundary conditions 


= ±1) = D$ q (z = ±1) = 0 


(7) 


for all q — 0, ±1,.... In the case q = 0 this is at variance with other authors 


13|, [14], [L6j who prefer 


D$ q=0 (z = ±1) = D 2 $ q=0 (z = ±1) = 0 


( 8 ) 


We think the latter condition is inconsistent in the nonlinear case where it 
would lead to a mixing of functions with different boundary conditions. As is 
discussed, e.g. in |20[] the boundary conditions (^) implies that the mass flux 
is kept constant whereas the case (§) refers to an experiment with constant 
pressure gradient averaged over the channel width. 

For the normal form transformation to follow we need the linearized part 
of (|) in diagonal form, if possible. Therefore we expand the functions in 
terms of the eigenfunctions Q q:l/ of the Orr-Sommerfeld equation, henceforth 
abbreviated by OSF 


^ d {x,z,t) - y Tj Q ,Jt)<& q!l ,(z)exp(iqax) , 


(9) 




where v = 1 , 2 ,... numbers the different eigenfunctions for given wave number 
q = 0, ±1, ±2, ±3,.... 

We solve the Orr-Sommerfeld equation numerically with standard proce¬ 
dures proposed in |17|, |2l[ using Chebyshev polynomials as basis functions. We 
took 50 polynomials both for the symmetric and the antisymmetric OSF. As 
a check of our IMSL procedure we reproduced the eigenvalues given in |I7| 
for the case a = 1, R = 10000. 

The case q = 0 can be solved analytically. Because of our unified boundary 
conditions ([7D we now get only antisymmetric functions, namely 

®q,v{z) — B v [s\n(k v z) - zsm(k„)] ; \ v = -k 2 u /R ; k u = tan(£y) . (10) 

In the usual case (^|) both symmetric and antisymmetric functions show up. 
However, one applies a kind of superselection rule suggested by numerical 
experience 


«VM = (-l)’ +1 4»,„(-i) , 

which excludes symmetric functions in the case q = 0. 
functions differ from the traditional ones 


( 11 ) 

Our antisymmetric 


z ) = C v [sm(vnz) + vnz(-l) u+1 ] ; A 0i! , = -u 2 n 2 /R . (12) 


4 





















Clearly, this does not affect the linear stability analysis because the critical 
amplitude is in the space with q = ±1. On the other hand, the different 
functions $ 0 ,^ should have an effect upon the nonlinear stability properties. 

Later on we measure the strength of the disturbance field in terms of the 
mean energy 


a 


Eh) = -u 2 + w 2 := — / dz / dx(u 2 + w 2 ) = 
2 87 r J -1 Jo 


— 2 E Vq,v{t)rjq'y{t)Mq 1/j qi l ,i , 
q,v,q' , 1 / 

which involves the metric components of the eigenfunctions as 


(13) 


Mq^q'v' — - 


-qq'oi 2 J dz$q !V $q>y + J dz 

where 6 denotes the Kronecker symbol. 


1 , d$q tV d$qiy 


dz dz 


bq'-q 1 ( 14 ) 


3 Galerkin representation 

In order to derive the Galerkin system for (£|) we need both the right and left 
eigenfunctions of the Orr-Sommerfeld equation. It is convenient to introduce 
the abbreviations j = (q,v), k = (<?', z/) and l = (r/', //'), and the index 
projection [j] = q, [k] = q' and [l] = q" and furthermore 

$*(*) = D$ k (z) ; G k (z) = (D 2 - [k} 2 a 2 )$ k (z) ; G k (z) = DG k (z) . (15) 

The Chebyshev expansion coefficients of $*,( 2 ), ®k(z), G k (z), G k (z ) are de¬ 
noted by a k , a k , g„, g k , respectively. By comparing the coefficients of the 
independent functions exp (iqax), we obtain 

E^r ~ X jVj](D 2 - q 2 a 2 )Qj(z) = iaJ2vkVi[^ki(z) - T k i(z)]8 [k]+[im (16) 
v 0,1 k,l 


where T ki (z) = [k\Q k (z)Gi(z) and T ki (z) = [l\<b k (z)Gi(z). The corresponding 
coefficients v kl and v kl are found according to the product rule given in [I7J 


v k J = 


M 

2 Cri 


OO 

E 

n'=—oo 


-k ~l ~kl 
a n-n'9n ' ! V n = 


J1 

2C77. 


OO 

■ E 

n'=— 00 


~k —l 

dn—n'9n' i 


n = 0,1,... (17) 


with a n = c\ n \a\ n \ and c n = 1 + S n:0 - Furthermore, the operator ( D 2 — q 2 a 2 ) is 
represented by the square matrix T> q as 

(D 2 - q 2 a 2 )T n (z) = V DJ,„T„,( Z ) , (18) 

n’ 
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where T n {z) denotes the Chebyshev polynomials of order n — 0,1,.... If q ^ 0, 
then T> g is nonsingular. Its elements can be inferred from formula (A5) given 

m 


m 


Comparing the coefficients of T n (z) in (|16|) we can write 

YJrfa - X i r U\Y, V nn' a h' = “ v* l }5 [k]+[l] . 


bl 


(19) 


k,l 


In the case q ^0 we multiply (|H1|) from the left by ( T> q )~ l and then by the 1- 
row matrix ti formed by the Chebyshev expansion coefficients b^n — 0,1,... of 
the left eigenfunction <E>h( z ). This gives rise to the nonlinear Galerkin system 
we are looking for 

-nVj = X jVj + J2 N mVkVi , (20) 

k,l 


with the constant coefficients 


Nj\ki = ia 


. En KEn'CP 


9 ^ — 1 \nM 

nn 1 Id n' 


V, 


kh 


dfci+ra.t?]; 


Em biiClri 

q + 0 ; j = (q, v) ; k = (s ', u ') ; l = (?", u ") ■ 


( 21 ) 


The denominator is obtained from the left and right eigenfunctions as delivered 
by the IMSL routine. 

In the case q — 0 the left and the right eigenfunctions <&o tV (z), given through 
(|T0|), are the same. We denote the corresponding expansion coefficients by d J n 
with j = (0, v). The Nj\ki are now given by 


N w = ™ : " l w ; j = ( 0,1/). (22) 

/ Jn dndn 

Of course, for all j only the symmetric part (Nj + Nj\ik) /2 enters (|2U|). By 
the reality of the Navier-Stokes equations we have the property 

N m= N w> < 23 ) 

where the star denotes complex conjugation and (q,u) = (—q,u). As an im¬ 
plication the amplitudes rjj obey the relation 


m = V*j • (24) 

By the Galerkin form ( |20j ) we tacitly assumed that the linear problem can 
be fully diagonalized. A degeneracy can be avoided by slightly shifting the 
Reynolds number R. 
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4 Interval normal form 


By means of a normal form transformation r/ —> £ we now decouple the 
Galerkin system © into a low-dimensiomal dominant system and a slaved 
subspace. The nonlinear, near identity transformation is written in multi¬ 
index notation as 

Vj = + £ A A m )^ m ; A A m ) e c ; i = 1,2,..., n , (25) 

meM 

where is a monomial of order \m\ = rrii + m 2 + ... + m n 

and M is the set of non-negative integer vectors 

M := {m\ m G Nq ; \m\ > 2} . (26) 

The transformed version of the Galerkin system then has the general form 

(l = Vf + £ . (27) 

meM 

which we try to simplify in the spirit of the normal form transformation as far 
as possible. In practical situations like our present case it is important to avoid 
not only zero denominators but also near resonances, otherwise the definition 
domain of the transformation (^) could become unacceptably small. I 11 view 
of this we introduce an interval parameter e and define the generalized resonant 
set 7 Zj(e). It determines the set of monomials to be kept in (|27|) as follows 

1Zj(e ) := {m \ dj(m ) < e, m G M} . (28) 

Here dj(m) is the absolute magnitude of a denominator of the normal form 
scheme 


Dj(m ) := (m, A) — A j ; (m, A) = miAi + rri 2\2 + ... + m n X n ; m G M . (29) 

The coefficients Aj(m ) and Bj{m ) of (^3]) and of the normal form © , respec¬ 
tively, are connected by recurrence relations of the following structure 

1 


A A m ) = 


Dj (m) 


[Bj(m) + P(A k (m'), Bi(m "))], 


(30) 


where P is a polynomial which contains only coefficients A k (m!) and Bi{m") of 
lower order rn 1 \, |m /7 | < \m\. If m G 7 Zj(e) then dj(m) = \Dj(m)\ < e and one 
sets the square bracket in ( |30| ) equal to zero with Bj{m) = —P(A k (m , ) } Bi(m ")). 
Otherwise Bj{m) is set equal to zero. In Appendix A we state the explicit re¬ 
lations for \m\ — 2 and \m\ = 3; the general recurrence relations can be found 
in 0. 

In qualitative normal form studies, see e.g. [f|, 0. fh| the interval pa¬ 
rameter e is set equal to zero. With e ^ 0 we now get the more complicated 
normal form 

•C, A,A, - £ B,(m)C' • (31) 

m&'R.j (e) 
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with more nonlinear monomials as compared with the pure case e = 0 . 

Nevertheless, with a suitable choice of e the interval normal form ([H]) ex¬ 
hibits a useful decoupling property. To reveal this, we partition the set of 
eigenvalues A j into groups characterized by their real parts. To this end the 
left hand complex plane of eigenvalues is divided into stripes which are paral¬ 
lel to the imaginary axis and of different width according to some convenient 
choice. The first stripe contains the least damped eigenvalues, i.e. with small¬ 
est modulus of the real part of A j. The second stripe is empty. The third 
stripe contains the second eigenvalue group. The fourth interval is empty, and 
so on. For illustration see fig. 1. 


Here insert fig. 1 


We label the groups by the index a and the members of one group by the 
index r. As should be noted, we have now three different equivalent labellings 
for the OSF: i) j = (q,v) where q = 0,±1,±2,... refers to the different 
streamwise wave numbers and v = 1, 2 ,... numbers the different states at a 
given q\ ii) j = 1,2,... which numbers the different OSF ordered e.g. by 
decreasing real parts 9?(Aj); iii) j = (ur) which is the labelling regarding 
decoupling, specifically 

n 

A j = Par + hu CTT ; a = 1,..., n ; r = 1,..., E a ; j = 1,..., ^ E a = n . (32) 

< 7 =1 

Note that the real parts p aT < 0 in this study. By convention, r = 1 and 
t = E a denote the maximum and minimum of the real parts for a given group 
a, respectively. The same labelling is carried over to rn 3 —»• rn aT and —*• 

Furthermore, we define the subspaces which correspond to the eigenvalues 
contained in one group: 


Sa := (£crl,£<T2, ■ ■ ■ ,£crE a ) ■ (33) 

We now define the width W a of a group by 

Wi = \piE 1 \ 

W a = | PaE a | - |Perl I for cr > 2 (34) 

and the distance A a between consecutive stripes containing eigenvalues by 

A<r IPct_i_i.iI IPctEctI • (35) 

Furthermore A a denotes an E a x E a matrix and P a a polynomial which starts 
with quadratic terms or is identically zero. As it turns out, the interval nor¬ 
mal form decouples substantially provided the distances A a between different 


groups are sufficiently large and the width W a of a group is sufficiently small. 
The precise theorem reads 

Decoupling theorem: If the interval parameter e allows for the properties 


> e - |p 
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and 


= A n 


^( 7 - 1 , 


T = Ai 


+ Pa 


for a = 1 , 2 ,... , 


(36) 

for (7 = 2,3,... , 


(37) 

llowing structure 



-i + A(=i) ; 


(38) 

“1, • ■ ■ ) “IT—1 ) j (7 = 2,.. 

.,h . 

(39) 


The proof of this theorem is given in Appendix B. As should be remarked, 
in practical application it is advantageous to consider also the imaginary parts 
of the eigenvalues, which leads to a smaller resonance set 7 Zj(e) and thus 
simplifies the normal form, in general. 

The structure of (^), (^) implies that the equations of the higher sub¬ 
spaces are subsequently linear with time dependent coefficients, if the solutions 
of the preceding subspaces are known. Moreover, if the subsystem Eh has a sta¬ 
ble fixed point at zero, then the fixed point is globally attractive for arbitrary 
initial values of the higher subspaces S^t = 0 ) with a = 2,... ,n provided 
the initial values of the subspace Si are in the basin of attraction of its fixed 
point. This is exemplarily proved in Appendix C by means of a Lyapunov- 
type function. The theorem is formal as long as the existence domain of the 


transformation (|25| ) is not established. If all eigenvalues are damped, then a 
finite convergence radius of (|25|) and thus a finite domain with Jacobian J ^ 0 
can be expected; in || this was proved for a special choice of e ^ 0 . 

In the subcritical region with negative real parts, the condition (^) can 
always be fulfilled, provided the interval parameter e is chosen small enough 
as 0 < e < |pn|. In this case (|37|) is unnecessary, because there is no condition 
on the mutual distances between different groups so that we can always choose 
W a = 0. This situation was considered in |J where the system Si is shown 
to contain only quadratic nonlinearities. However, with such a choice the 
parameter e would become arbitrarily small in the limit R —> R c , with the 
implication that the definition domain of (^5|) may shrink to zero. 

In our application we choose e > |pn|. The problem is now hireling a 
suitable partitioning into groups which is compatible with the conditions (j36l), 
Moreover, it is desirable to maintain a feasible form of the dominant 
subsystem Si. With the given finite cut-off of 22 eigenfunctions we have no 
problems finding a convenient partitioning. 
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5 Application to Poiseuille flow 

We choose the interval parameter e as 

e = max{ | X q=0 ,u=i \, 21 \ q=ltV=1 ) |} (40) 

e < |3 : (A g= i i j /= i) | , 

where 3ft and 9 denote the real and imaginary part, respectively. This choice 
provides us with a three-dimensional dominant subsystem in the Reynolds 
number interval 4320 < R < R c . The partitioning of eigenvalues into groups 
and with it the structure of the normal form is the same in this R-interval. The 
situation is quantitatively shown in fig. 1 for R = 4320. It is compatible with 
the conditions ([I6|) and (0) for the distances A CT and widths W a of the stripes. 
As should be noticed, our chosen interval parameter e increases with decreasing 
Reynolds number R according to the R dependence of the eigenvalues. 

The first stripe contains one real and the two complex conjugate eigenval¬ 
ues which become critical at R = R c . For convenience we denote the three 
dominant variables and eigenvalues by (£i,£ 2 , £ 3 ) and (Ai, A 2 , A 3 ), respectively. 
Because Ai = Ag and A 3 = A 3 we have the properties £1 = ££ and £3 = £ 3 - 


In view of the normal form (|3T|) we introduce the abbreviations 

r 2 = £16 > 0 (41) 

[A] 

dji 7 " 2 ) = B j ( n /B + n Afc + fij)r 2n ; N = -- ; J = 1,2,3 (42) 

n =1 2 1 pul 

[A] _ , | x i 

h(r 2 ) = ^ 3(^1 + nix 2 )r 2n ; N = ' " ; (43) 

n= 1 21 Pii| 

a*, = 2W-13 ; k — 1,2 (44) 

where [N] and [N] denote the integer parts of N and N, respectively. With 
this, in the subspace Si the normal form ([JI]) has the following structure 

£1 = [Ai + 016 +gi(r 2 )]£i 

£2 = [A 2 + a 2 £ 3 + ^(r 2 )]^ (45) 

£3 = A 3 £ 3 + h{r 2 ) + p 3 (r 2 )£ 3 


Obviously, the above system is decoupled from the remaining variables with 
k > 4. 

In the following we take the transformation (|25|) and with it the normal 
form (|3T|) up to cubic terms. With the aid of polar coordinates £1 = rexp(z<p) 
and by setting £ 3 = £, A 3 = A, pn = p and Uu = oj the transformed dominant 
system is further reduced to effectively two dimensions 

r = r[p + 9?(cii)£ + br 2 ] ; r>0; (46) 
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£ = ^ + lir 2 + 12T 2 ^ ; 
<fi = lu + /3r 2 . 


(47) 

(48) 


As illustrated in fig. 2, we obtain a stable fixed point at r = £ = 0 and a 
further fixed point A* at (r c , £ c ), which corresponds to a Hopf limit cycle |23 
with angular frequency given by (S3). 


Near the critical point R = R c we have £ c = — 7 ir 2 /A which renormalizes 
([46]) to r = r[p + 6i(i? c )r 2 ] with the Landau-type parameter at R = R c 


tl = 6- 7l * (ai > . (49) 

A 

For higher order dominant systems there enter further renormalizing sum¬ 
mands into bi- 

From table 1 it is seen, that the Landau parameter 6^ of our normal form 
method is almost independent of the choice of the interval parameter e, as it 
should be. 


Table 1 

Dependence of the Landau parameter bi on the choice of the interval parameter 
e at R = 5770 with E\ denoting the dimension of the corresponding dominant 
system. The main results of the paper are for E\ = 3. 


e 

A g=0)Z ,=2 

V 

kQ 

II 

o 

II 

2|K(A 9=1>i/=1 )| 

b L 

39.163 

39.168 

39.174 

E\ 

5 

3 

2 


As can be inferred from fig. 2, the basin of attraction of the zero fixed 
point is limited by the curve r — 0. If the starting point lies inside this basin 
of attraction, both the variables of the dominant and the slaved system go to 
zero. 


Here insert fig. 2 


Let us examine the slaved system at the non-zero fixed point (r, £) = 
( r e ; £c)- To this end we write down in the following the full normal form of 
our Galerkin system at R = 4320 with the first two equations referring to the 
dominant system (|46|), (|47|) 

f = r(p + 0.02942£ + 46.64r 2 ) ; r>0 
£ = £(A - 2.568r 2 ) + 0.6561r 2 
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(50) 

(51) 











£4 

= £ 4 (A 4 - 58.43r 2 ) - 75.72£r 2 

(52) 

4 

= £ 5 (A 5 - 110.19r 2 ) - 65.94£ 4 r 2 

(53) 


= f 6 (A 6 + 80.37r 2 ) 

(54) 

r 1 

= r 1 (pi — 0.06202£ + 12.39r 2 ) ; r 4 > 0 

(55) 

6 

= £ 9 (A 9 + 388.6r 2 ) 

(56) 

r 2 

= r 2 (p 2 - 0.08132£ + 83.36r 2 ) ; r 2 > 0 

(57) 

r 3 

= r 3 (p 3 - 0.09851£ - 206.2r 2 ) ; r 3 > 0 

(58) 

£14 

= ^ 4 (Ai 4 + 168.2r 2 ) 

(59) 

£15 

= £is(Ai 5 - 155.95r 2 ) + 2561£ 4 r 2 

(60) 

Cl6 

= U(Ai6 - 61.77r 2 ) - 86210^2 

(61) 

£l7 

= 6t(Ai 7 + 14.79r 2 ) + 15619£ 5 r- 3 2 + 63.38£i 4 r 2 

(62) 

£l8 

= 6s(Ai 8 + 0.6672r 2 ) + 19360£ 44 r 2 + 5671£ 9 r 2 

(63) 

£l9 

= 6 9 (Ai 9 - 1.430r 2 ) + 1280^i5r 3 

(64) 

£20 

= 6o(A 20 - 1.772r 2 ) - 4.482£ 18 r? - 238.51£ 17 r 2 2 

(65) 

^21 

= ^ 2 i(A 2 i — 1.095r 2 ) — 0.7456£i 9 r 2 

(66) 

£22 

= 6 2 (A 22 - 0.3739r 2 ) - 172.1£ 19 r 2 , 

(67) 

magnitudes p = 9?(Ai), A = A 3 , A 4 ,..., A 22 of the linear problem are 


listed in Appendix D. As in the case of the dominant system, polar coordinates 
are used for the complex conjugate amplitudes with r i, r 2 and r 3 denoting the 
moduli of the modes (q = ±1, v = 2), (q = ±2, v = 1) and (q = ±3, v = 1), 
respectively. The equations for the corresponding polar angles are omitted. 

have nonzero fixed points 


The equations fl52|) and (|53|) have nonzero fixed points and £ 5 , respectively, 
whereas each of the remaining equations has fixed point zero. Thus the overall 
nonzero fixed point at R = 4320 is given by 

(r c = 0.008231, = 0.009170, = -0.002647, 

C 5 = 0.0003380, a = 0, • • •, £22 = 0 ) . (68) 

When the variables r and £ are substituted by r c and £ c in the slaved equations 
©-©. then it turns out that the separate fixed point of the slaved system 
is stable. Therefore, at the critical point of the dominant subsystem there is 
no runaway of the slaved variables. We can thus define a meaningful critical 
energy by mapping the coordinates of the nonzero fixed point back into the 
original phase space of the amplitudes p. With the aid of (|13|) we obtain 
E = 0.26 * 10~ 3 at R = 4320 which has to be compared with the unperturbed 
energy Eq = 8/15 of the basic flow. The relative critical energy E / E 3 has the 
order of magnitude established in the literature, see e.g. 
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Regarding the fixed point in the original phase space, we state the maximal 
amplitude for a given wave number q 


r/"“o = 9.176 * 10 


-3 


r/™ff = 5.726 * 10 


-3 
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|?/™ = a 2 T | = 2.630 * 10~ 8 ; \r]™*\ = 3.577 * 10" 9 . (69) 

As a final remark of this section, it is straightforward to extend the given 
formalism to Reynolds numbers below R = 4320 by choosing a larger interval 
parameter e. This amounts to adding further nonlinear terms to the normal 
form so as to preserve the same structure in the larger Reynolds number in¬ 
terval. 


6 Comparison with Landau method 


Our main results are represented by fig. 3. We adopt the same normalization 


$,=!,«/=! (* = 0 ) = 1 


as used in the Landau method [O, [12 


(70) 


Here insert fig. 3 


As is seen both the critical curve of the Landau method and of our normal 
form approach are consistent with the neutral curve of the Galerkin system, 
which is approximated by directly integrating the Galerkin system. The start¬ 
ing points were choosen close to the critical amplitude vector r) s , which is the 
image of the fixed point F* of the normal form space. This result also indicates 
that, in view of the Galerkin system choosen, our function space of 22 Orr- 
Sommerfeld eigenfunctions is sufficiently large. However the Galerkin space is 
to small to reproduce the Landau method results near the critical point for 
5760 < R < R c . We find b L (R c ) = 39.17, which is by 36% larger than the 
corresponding value, 28.88, obtained by the Landau method |12[]; correspond¬ 
ingly at the critical Reynolds number our critical amplitude is by 16% smaller, 
see fig. 3 and fig. 4. As should be remarked the critical curve of the Landau 
method |T^| is established in the smaller interval 5300 < R < R c . 

For comparison with the results reported elsewhere [|l^] we give the angular 
frequency of the nonlinear wave at the fixed point F* for several Reynolds 
numbers. As table 2 shows, the numerical values agree within 1%. 


Table 2 

Comparison of the angular frequency of the nonlinear wave with the Landau 
method. 


R 

5200 

5400 

5500 

5700 

5770 

F 

Landau 

0.2771 

0.2743 

0.2744 

0.2725 

0.2730 

0.2717 

0.2704 

0.2699 

0.2694 

0.2693 
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In fig. 4 we show results for the parameter 6^ as a function of the number 
N of modes with q — 0 considered. The new boundary conditions, curve (2), 
and the usual ones, curve (1), give practically the same value b^. 


Here insert fig. 4 


As to the definition domain of our normal form transformation we cannot 
safely rely on the criterion of the Jacobian. In the given case the Jacobian is 
larger 1 well beyond the critical amplitude. 


7 Conclusions 

A generalized normal form scheme for dissipative flows was presented which 
allows to control the minimal magnitude of the normal form denominators, 
and at the same time keeps a useful decoupling property of the normal form. 
The decoupling property was rigorously proved and constitutes a new result. 
The feasibility of the method was examined for the well studied example of 
plane Poiseuille flow in the subcritical interval 4320 < R < R c of the Reynolds 
number R. The critical energy E C (R) or equivalently the critical amplitude 
rj c (R), at which the laminar flow becomes unstable, agrees well with the results 
obtained by other methods |T2|, except close to the critical point where our 


critical amplitude is about 16% smaller than the corresponding value of the 
Landau method |nj. The latter is established in the interval 5300 < R < R I: . 
whereas the present normal form gives reliable results down to R = 4320,0 see 
fig- 3. 

The function space, in our case the Orr-Sommerfeld eigenfunctions, is rig¬ 
orously determined in principle for every parameter point R. For comparison, 
both the Landau and the center manifold method rely on a critical parameter 
point R c of the linearized problem around which one expands in powers of 
y/R c — R, for instance. The normal form method, on the other hand, requires 
the existence of a stable stationary solution only which provides a fixed point 
in phase space. Therefore the present method should be useful to study e. g. 
pipe flow or plane Couette flow for which the linearized system lacks a critical 
point. 

As a minor detail we studied boundary conditions which differ from the 
usual ones for the modes with q = 0. Our new boundary conditions constitute 
a unified function space for all modes. However it turns out that both boundary 
conditions give practically the same results, see fig.4. 

The quantitative validation of an existence domain is a notoriously difficult 
problem in all three discussed methods. We found our normal form results be¬ 
ing consistent with results from the direct integration of the Galerkin system. 
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At the critical curve r c (R) we found the Jacobian of the normal form trans¬ 
formation being always positive and well away from zero. However, since our 
normal form transformation neglects higher than cubic terms, the Jacobian 
criterion is not safe here. 

The numerical efforts reside to about 90% in the normal form method, 
whereas the remaining computing time lies mainly in establishing the Orr- 
Sommerfcld eigenfunctions and the Galerkin system. For one parameter value 
R the computing time was about six minutes on a RISC 6000 work station. 

The parameter interval 4320 < R < R c considered was determined by 
the choice of the interval parameter e together with the requirement that the 
normal form (|3l|) keeps the same structure in the whole parameter interval. It 
is in principle possible to go deeper into the subcritical range by choosing a 
larger e which would result in a higher dimensional dominant subsystem. 


Appendix A: Recurrence relations 


We insert the normal form transformation (25) into the Galerkin system (20) 
and consider the normalized dynamical system (|27| ) as an ansatz. The left 
hand side of (|20|) then reads 


Vj = ^ [(m, X)Aj(m) + Bj(m )]f m 

mGM 


+ - m'+ , (71) 

m,m'£M l 


where fx k denotes the unit vector with (// fc )j = 5 ik . We use the convention that 
Aj(m ), Bj(m ) = 0 if m does not belong to the basic set M defined in 
The right hand side of ( |20|) becomes 


Aj£j + Xj ^2 A'( m )£ m + Nj\ki 

m£M k,l 




H [M m - hO 

mGM 


+A k (m — Hi)\£ m + J2 A fc (m - m')Ai(m')^ m \ ■ (72) 

m,m'£M J 

By equating the two sides, the linear terms cancel out. The coefficients of the 
second order terms with \m\ — 2 read 

Aj{m) = | N j\ki8( m i hfc + hO if rn 7 Zj(e) ; 

^jv 71 ) k,l 

Bj(m ) = -D j (m)A j (m) + N j\kiK m , Vk + Hi) if me fcj(e) , (73) 

k,l 
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where 5(m,m') = 1 for m = ml and 5(m,m') = 0 for m 7 ^ ml. Note that 
Aj(rn) can be chosen arbitrarily if m G lZj(e). In the case \m\ = 3 we get 


Aj(m) 


1 

Dj (m) 


2 J2 N j\kiA k (m - in) 

k,l 


B k {m’) [m k - m' k + 1 ]Aj(m - m + Pk) 

k=l m'£ 7£fc(e) 


if m IZj (e) ; 


Bj(m) = -Dj{m)Aj{m) + 2 J2 Nj\kiA k (m - pi) 

k,l 


n 

B k (m')[m k - m' k + l]Aj(m - m'+ fi k ) if melZj(e). (74) 

k =1 m' ^'JZk(e) 

Once more Aj(m) can be chosen arbitrarily if m G IZj(e). 

For convenience we set the resonance coefficients Aj{m) = 0. 


Appendix B: Proof of decoupling theorem 

The forms (^) and (|39|) of the normalized Galerkin system imply that the 
dynamics of the subspace is independent of the higher subspaces Ev with 
cr' > cr and depends only linearly on its own variables if a > 2. We prove 
this by the method of contradiction under the conditions m-m Let us 
consider the system S CT , and assume that i) in the normal form (|3l|) there is 
a coefficient B aT {m) 7 ^ 0 with m a ' T ' 7 ^ 0 and a' > a or that ii) there exists at 
least one coefficient B ar (m) 7 ^ 0 with YJl?=\ > 2 for a > 2. In case i), 
because \m\ > 2 always, and by ([36]) the estimates are 

dor(m) A |(m)) | k |p cr / r / pw Par \ 

| Pa'r' Pcrr| 4“ |Pll| ^ T |Pll| ^ i (75) 

which tells that m ^ 77j(e) with j = (cr, r) and thus B aT {m ) = 0 in contradic¬ 
tion to the assumption i). 

In the case ii) we get at first 

(7 E f 

d aT {m ) > |9?(A rr (m))| > I J2 £ r'Pc'r’ - par] • (76) 


a 

I EE r kTl(j , T l Pcr'r' | ^ | P (TT (77) 

a'= lr'=l 

then by diminishing the left hand side under the assumption ii) we can write 

a E a r 

2\Pal\ < | Par' \ + \Par»\ < EE ^a'r' Pa'r' | ^ \Par\ — I/LtEctI (78) 

cr '=1 r '=1 
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which implies \p a \\ < W a and contradicts (|37|). In the case complementary to 
(|77|) we have 


d ryT {lfTl) ^ A ' y ' Tfl a ' T i\Pa’r' | |/7rr| £ 2|p CT i| \p&E a \ • (79) 

cr'=l t'—\ 

Making use of (^7|) in the form \p a i\ > W a + e we continue (|79|) as 

d aT (m) > \p al \ + W a + e - \paE a \ = e (80) 

which means that case ii) is not contained in the resonant set 7 Zj(e) with 
j = (cr, r). This completes the proof. 


Appendix C: Proof of slaving 

The dominant subspace Si is independent of the higher subspaces S CT , a > 2. 
Since the linear part of S] is stable, the flow S x possesses a finite basin of 
attraction BA by the theorem of Hartmann-Grobman ||. We claim that, 
if the initial points Si (f = 0) G BA , then the fixed point £ = 0 is stable 
independently of the initial values S a {t = 0) with a > 2. We demonstrate 
this in an exemplary way for a three-dimensional system. The general proof 
merely needs more efforts in writing, but is otherwise fully analogous. 

Let us assume that the first subspaces Si(t),..., S (r _i(t) are stable, and 
that the next subspace S CT is three-dimensional with one real variable x and 
two complex conjugate variables rexp(z^), rexp(— itp) with r > 0. By ( |5U| ) 
the corresponding dynamical system reads 

x = [A^i + 7 i (t)]x + 7 2 (t)r + 0 i(t) ; (81) 

r = 7 3 (t)x + [p a2 + 74 (t)}r + 6 2 (t) , (82) 

where both p a2 = ^(Ao- 2 ) and Ao-i are negative, and the real coefficients 7 (t), 
&(t ) are determined by the variables of the lower subspaces. We need not know 
the phase ip = yit). By the inductive assumption, the four functions 71 to 74 
and 0 i /2 become arbitrarily small for t —> 00. Since the above system is linear 
with regular time dependent coefficients, the solutions x(t) and r(t) certainly 
are bounded within any finite time interval. With A = maxIAo-i, p a2 } < 0 the 
time derivative of the function £ = x 2 + r 2 obeys the relation 

C/2 < X[x 2 + r 2 ] + 7 i £ 2 + y 4 r 2 + [72 + 7 3 \xr + Q x x + 0 2 r . (83) 

Applying an orthogonal transformation (x,r) —> (yi,y 2 ) we can get rid of the 
mixed term xr. The transformed expression reads 

C/2 < A ’[yl + y^\ + d\y\ + id 2 y 2 , (84) 
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with X' — X + 5(t) where S(t) and $ 1/2 are arbitrarily small magnitudes if t is 
sufficiently large. As a consequence, if we choose points with 


tel> W ; i = 1 ’ 2 ’ (85) 

then £. < 0 and all initial points end asymptotically in the complementary, 
arbitrarily small domain T> 0 with y 3 \ < |2i? J -/A , |. As a final step we invoke the 
theorem of Hartman-Grobman |5J which states that our transformed system 
(J38|) , ( |39D is homeomorph to the linearized problem in a finite neighbourhood of 
the origin, provided all eigenvalues have a negative real part. If t is sufficiently 
large, then T >0 is a subset of the validity domain of the Hartman-Grobman 
theorem. Thus, if t —> 00 then (x,r) —> (0,0). 

Appendix D: Table of eigenvalues 

Table 2 

Real and imaginary part 3?(A), A(A) of the eigenvalues X 3 with j = (q,u) for 
R = 4320 and R = 5000 at a = a c . The ordering is given according to the 
real parts at R — 4320. The different eigenvalues at given wave number q 
are labeled by v — 1,2,.... S and A denote symmetric and antisymmetric 
eigenfunctions, respectively. 
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No. 

Q 

V 

S/A 

Sft(A) 

R = 4320 

A(A) 

R = 4320 

r 9?(A) 

R = 5000 

A(A) 

R = 5000 

1 

- 1 

1 

S 

-0.0034297871 

0.28323209 

-0.0015441660 

0.27621304 

2 

1 

1 

S 

-0.0034297871 

-0.28323209 

-0.0015441660 

-0.27621304 

3 

0 

1 

A 

-0.0046737834 

0 

-0.0040381470 

0 

4 

0 

2 

A 

-0.013814702 

0 

-0.011935903 

0 

5 

0 

3 

A 

-0.027523116 

0 

-0.023779971 

0 

6 

0 

4 

A 

-0.045800419 

0 

-0.039571561 

0 

7 

- 1 

2 

S 

-0.053931898 

0.96621013 

-0.050157877 

0.97004199 

8 

- 1 

2 

S 

-0.053931898 

-0.96621013 

-0.050157877 

-0.97004199 

9 

0 

5 

A 

-0.068646855 

0 

-0.059310882 

0 

10 

- 2 

1 

A 

-0.075518479 

1.9643698 

-0.070275192 

1.9697702 

11 

2 

1 

A 

-0.075518479 

-1.9643698 

-0.070275192 

-1.9697702 

12 

- 3 

1 

S 

-0.090121159 

2.9671465 

-0.084056136 

2.9738632 

13 

3 

1 

S 

-0.090121159 

-2.9671465 

-0.084056136 

-2.9738632 

14 

0 

6 

A 

-0.096062496 

0 

-0.082997997 

0 

15 

0 

7 

A 

-0.12804737 

0 

-0.11063292 

0 

16 

0 

8 

A 

-0.16460149 

0 

-0.14221568 

0 

17 

0 

9 

A 

-0.20572486 

0 

-0.17774628 

0 

18 

0 

10 

A 

-0.25141749 

0 

-0.21722471 

0 

19 

0 

11 

A 

-0.30167938 

0 

-0.26065098 

0 

20 

0 

12 

A 

-0.35651052 

0 

-0.30802509 

0 

21 

0 

13 

A 

-0.41591092 

0 

-0.35934704 

0 

22 

0 

14 

A 

-0.47988059 

0 

-0.41461683 

0 
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Figure captions 


Fig. 1. Partitioning of the eigenvalues into groups for the de¬ 
coupling theorem at R = 4320. The eigenvalues are indicated by 
circles. The first stripe of width W\ defines the dominant system; 
it contains the two complex conjugate eigenvalues which become 
critical at R = R c and one real eigenvalue. All other stripes con¬ 
taining eigenvalues have zero width here. The distances are 
all larger e — \pn\ : where e is the interval parameter of the nor¬ 
mal form. Note that one easily could group the eigenvalues in a 
different way when the distances become too small. 

Fig. 2. Critical lines r = 0 and £ = 0 in £r-space as given 
by the transformed dominant system (QT|), ([47]) for R = 4320. 
The stable fixed point is at the origin (0,0), the unstable one is 
denoted by F*. The arrows indicate the direction (r, £) of the 
flow. For r > 0.0427, which is out of scale, the arrows along the 
curve £ = 0 change sign. 

Fig. 3. Critical amplitude r] c (solid curve) by the normal form 
method as a function of the Reynolds number R , and stability 
boundaries found by direct integration of the Galerkin system 
(dashed curves); the up and down arrows indicate the existence 
of growing modes and decay of all modes, respectively. The criti¬ 
cal amplitudes of the normal form (1) and of the Landau method 
(2) are compared in the insert. 

Fig. 4. Dependence of the renormalized parameter bi on the 
number N of the eigenfunctions with q = 0 for R = 5770. Curve 
(1) and curve (2) refer to the usual and our boundary conditions, 
respectively. In both cases the Orr-Sommerfeld basis functions 
are normalized to the same energy. 
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